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Abstract —Nonlinear and dispersive transmission impairments 
in coherent fiber-optic communication systems are often com¬ 
pensated by reverting the nonlinear Schrodinger equation, which 
describes the evolution of the signal in the link, numerically. This 
technique is known as digital backpropagation. Typical digital 
backpropagation algorithms are based on split-step Fourier 
methods in which the signal has to be discretized in time and 
space. The need to discretize in both time and space however 
makes the real-time implementation of digital backpropagation 
a challenging problem. In this paper, a new fast algorithm for 
digital backpropagation based on nonlinear Fourier transforms is 
presented. Aiming at a proof of concept, the main emphasis will 
be put on fibers with normal dispersion in order to avoid the issue 
of solitonic components in the signal. However, it is demonstrated 
that the algorithm also works for anomalous dispersion if the 
signal power is low enough. Since the spatial evolution of a 
signal governed by the nonlinear Schrodinger equation can be 
reverted analytically in the nonlinear Fourier domain through 
simple phase-shifts, there is no need to discretize the spatial 
domain. The proposed algorithm requires only 0{D\og^ D) 
fioating point operations to backpropagate a signal given by D 
samples, independently of the fiber’s length, and is therefore 
highly promising for real-time implementations. The merits of 
this new approach are illustrated through numerical simulations. 

Index Terms —Optical Fiber, Digital Backpropagation, Nonlin¬ 
ear Fourier Transform, Nonlinear Schrodinger Equation 

I. Introduction 

The evolution of a signal E{x^t) e Cm optical fiber, where 
X ^ 0 denotes the position in the fiber and t g M the time, 
is well-described through the nonlinear Schrodinger equation 
(NSE). Through proper scaling and coordinate transforms, the 
NSE can be brought into its normalized form 

ij^ + ^+2K\E\^E = -iVE, K6{+1}. (1) 

The parameter n in Eq. 0 effectively determines whether the 
fiber dispersion is normal (—1) or anomalous (+1), while the 
parameter T determines the loss in the fiber. In the following, it 
will be assumed that the loss parameter is zero. There are two 
important cases in which a fiber-optic communication channel 
can be modeled under this assumption. When the fiber loss 
is mitigated through periodic amplification of the signal, the 
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average of the properly transformed signal is known to satisfy 
the NSE with zero loss m . Eurthermore, recently a distributed 
amplification scheme with an effectively unattenuated optical 
signal (quasi-lossless transmission directly described by the 
lossless NSE) has been demonstrated |i2l. If the loss parameter 
is zero, the NSE can be solved using nonlinear Fourier 
transforms (NFTs) ID The spatial evolution of the signal 
E{x^t) then reduces to a simple phase-shift in the nonlinear 
Fourier domain (NFD), similar to how linear convolutions 
reduce to phase-shifts in the conventional Eourier domain. The 
prospect of an optical communication scheme that inherently 
copes with the nonlinearity of the fiber has recently led to 
several investigations on how data can be transmitted in the 
NED instead of the conventional Eourier or time domains 

0, O, El, Q, il, a, Col, CD, with the original idea 

being due to Hasegawa and Nyu 1121 . Next to potential 
savings in computational complexity, it is anticipated that 
subchannels defined in the NED will not suffer from intra¬ 
channel interference, which is currently limiting the data rates 
achievable by wavelength-division multiplexing systems CD- 
The mathematics behind the NET are however quite in¬ 
volved, and despite recent progress in the implementation of 
fast forward and inverse NETs C3, ca, ca no integrated 
concept for a computationally efficient fiber-optic transmission 
system that operates in the NED seems to be available. In 
this paper, therefore the problem of digital backpropagation 
(DBP), i.e. recovering the fiber input E(f)^t) from the output 
E{xi^f) by solving is addressed instead using NETs 0. 
Both concepts are compared in Eig. Although DBF does 
not solve the issue of intra-channel interference in fiber-optic 
networks because it is usually not feasible to join the individ¬ 
ual subchannels of the physically separated users into a single 
super-channel in this case ini X.D], we note that there are 
other scenarios where this issue does not arise lEl. The advan¬ 
tage of digital backpropagation over information transmission 
in the NED is that it is not necessary to implement full forward 
and backward NETs. Together with recent advantages made 
in ca, ca, oa, this observation will enable us to perform 
digital backpropagation in the NED using only 0{D\og^ D) 
fioating point operations (flops), where D is the number of 
samples. This complexity estimate is independent of the length 
of the fiber because the spatial evolution of the signal will 
be carried out analytically. Conventional split-step Eourier 



Figure 1. Information transmission (top) vs DBF (bottom) in the NFD 


methods, in contrast, have a complexity of O[MD log D) 
flops, where M is the number of spatial steps (Tsl III.G]. 

The goal of this paper is to present a new, fast algorithm 
for digital backpropagation that operates in the NFD and to 
compare it with traditional split-step Fourier methods through 
numerical simulations. The impact of noise resulting from the 
use of distributed Raman ampliflcation will be of particular 
interest. The paper is structured as follows. In Sec. [n| the 
theory behind digital backpropagation in the NFD will be 
outlined, and the new, fast algorithm will be given and 
discussed. The simulation setup is described in Sec. |I^ while 
results are reported in Sec. |IV] Sec. |V| concludes the paper. 

II. Digital Backpropagation in the 
Nonlinear Fourier Domain 


In this section, first the theoretical and computational results 
that are required to perform digital backpropagation in the 
NFD are briefly recapitulated from CS). Afterwards, the fast 
algorithm is presented and its limitations are discussed. 


A. Theory for the Continuous-Time Case 


The Zakharov-Shabat scattering problem associated to any 
signal E{x^f) that vanishes sufficiently fast for t ^ +oo is 
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With and 02 denoting the components of 0, now define 
a(x,X) := lim 0i(x, A), (4) 

t^oo 

0(x,A) := lim (j) 2 {x^t^X)^ 

t^oo 

where A g C is a parameter. If E{x^ t) satisfies the NSE 0 
with zero loss parameter F = 0, the corresponding a{x, A) 
and (3{x,X) turn out to depend on x in very simple way: 

(t(x, A) = (t( 0 , A), /3(x, A) = ^ 0(0, A). (5) 


These functions are not the final form of the NFT, but for our 
needs it will be sufficient to stop here. 

The complicated spatial evolution of the signal E{x,t) thus 
indeed becomes trivial if it is transformed into a{x,X) and 
0(x,A). Based on this insight, ideal continuous-time digital 
backpropagation reduces to three basic steps in the NFD: 
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B. Discretization of the Continuous-Time Problem 

In order to obtain a numerical approximation of a{x,X) 
and 0(x, A) for any fixed x, choose a sufficiently large interval 
[Ti, T 2 ] inside which E{x, f) has already vanished sufficiently. 
Without loss of generality, one can assume that Ti = — 1 and 
T 2 = 0. In this interval, D rescaled samples 


E[x^ n] := eE ■ 


-l + ne-- 


^ , e n s 


are taken. With z := e ([^-([^ becomes 

1 z~^E[x,n] 


0 [x, n, z] :=z‘^ 


(f)[x,0, z] :=z 2 


-K,E[x,n] z 
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This leads to the following polynomial approximations: 

a{x,X) % a{x,z) '=T^fjo^ai{x)z~'' := 0 i[x,D, 2 :], 
l3{x,X) ^ b{x,z) ■■=TifjQbi{x)z~^ := (p 2 [x,D,z]. 


C. The Algorithm 

The three steps in the diagram ^ can be implemented with 
an overall complexity of 0{D\o^D) flops as follows. 

Step A : The discrete scattering problem to And the poly¬ 
nomials a{xi^z) and Mxi^z) from the known fiber output 
E{xi^t) through (|^-([^ is solved with only 0{D\og^ D) 
flops using (151 Alg. 1] (also see (141 Alg. 1]). 

Step B: The unknown polynomials a(0, and hift^z) are 
defined by the unknown fiber input F^( 0 ,t) through 
In this step, approximations d(0, z) and 6(0, z) of a(0, z) and 
6 ( 0 , 2 ;) are computed based on Eq. The left-hand in ^ 
suggests the choice d(0, z) := a{xi,z). For finding 6(0, t), let 
us denote the D-th root of unity by w := . The right- 

hand in ^ motivates us to And 6 ( 0 , 2 ;) = bi{0)z~'^ by 

solving the well-posed interpolation problem 

=e4i(log(^"-i)/(-2ie))%, (9) 

= b{xi,w^-^),ne 


with the fast Fourier transform, using only 0{D\ogD) flops. 

Step C: The inverse scattering problem of estimating 
E[0^ n] from d( 0 , z) and 6 ( 0 , z) by inverting (|^-([^ is solved 
using 0{D log^ D) flops as described in (161 IV]. 


D. Limitations 

It was already mentioned above that the algorithm works for 
fibers with normal dispersion. In that case, the sign n in the 
NSE 0 will be negative and E{x^t) is determined through the 
values that q;( 0. A) and 0(0, A) take on the real axis M 9 A (191 
p. 285]. If the dispersion is anomalous, the sign will be positive 
and E{x, f) is determined through the values that q;( 0. A) and 
0(0, A) take on the real axis and around the roots of q;( 0. A) 
in the complex upper half-plane ^(A) > 0 (191 IV.B]. Taking 
the coordinate transform 2 ; = into account, one sees 
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Figure 2. a) Simulation setup of coherent optical communication systems 
with DBF in the NFD, b) basic block functions of OFDM and Nyquist 
transceivers 


that the interpolation problem ([^ however enforces the phase- 
shift only for certain real values of A. Thus, the algorithm is 
unlikely to work if ct(0, A) has roots in the upper half-plane. 
The condition \E{0^t)\dt < f is sufficient to ensure that 
there are no such roots ||20l Th. 4.2]. The actual threshold 
where solitons start to emerge however is expected to be higher 
due to the randomly oscillating character of waveforms used 
in fiber-optic communications 12T1 . (221. 

III. Simulation Setup 

In this section, the simulation setup that was used to assess 
the performance of DBF in the NFD is presented. 

The transmission link was assumed to be lossless due to 
ideal Raman amplification, with an amplified spontaneous 
emission (ASE) noise density A^ase = ^LH/sKt. Here, T 
is the fiber loss, L is the transmission distance, hfs is the 
photon energy, fg is the optical frequency of the Raman 
pump providing the distributed gain, and Kt = 1.13 is 
the photon occupancy factor for Raman amplification of a 
fiber-optic communication system at room temperature. In 
the simulations, it was assumed that the long-haul fiber link 
consisted of 80-km spans of fiber {standard single-mode fiber 
(SSMF) in the anomalous case) with a loss of 0.2 dB/km, 
nonlinearity coefficient of 1.22 W/km, and a dispersion of 
+ 16 ps/nm/km (normal and anomalous dispersions). A photon 
occupancy factor of 4 was used for more realistic conditions. 
The ASE noise was added after each fiber span. 

The data was modulated using high spectral efficiency mod¬ 
ulation formats (QPSK and 64QAM) and either Nyquist pulses 
(i.e. sine’s) M or orthogonal frequency division multiplexing 
(OFDM) (71, O. The block diagram of the simulation 
setup and basic block functions of the OFDM and Nyquist 
transceivers are presented in Fig. HI In Fig. I^a), DBF in the 
NFD is performed at the receiver after coherent detection, syn¬ 
chronization, windowing and frequency offset compensation. 
For simplicity, both perfect synchronization and frequency 
offset compensation were assumed. The net data rates of 
the considered transmission systems were, after removing 7% 
overhead due to the. forward error correction (FEC), 100 Gb/s 
and 300 Gb/s for QPSK and 64QAM, respectively. For the 
OFDM system, the size of the inverse fast Fourier transform 
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Figure 3. a) Illustration of a burst mode transmission at the transmitter side, in 
which neighboring packets are separated by a guard time, b) received signal 
at the receiver, which is broadened by chromatic dispersion, c) windowing 
signal for processing in the proposed nonlinear compensation scheme. 


(IFFT) was 128 samples, and 112 subcarriers were filled with 
data using Gray coding. The remaining subcarriers were set 
to zero. The useful OFDM symbol duration was 2 ns and 
no cyclic prefix was used for the linear dispersion removal. 
An oversampling factor of 8 was adopted resulting in a total 
simulation bandwidth of '-448GHz. At the receiver side, all 
digital signal processing operations were performed with the 
same sampling rate. The receiver’s bandwidth is assumed to 
be unlimited in order to estimate the achievable gain offered 
by the proposed DBP algorithm. 

The data was transmitted in a burst mode where data 
packets where separated by a guard time. See Fig. for 
an illustration. The guard time was chosen longer than the 
memory AT = 27rB/32L induced by the fiber chromatic 
dispersion (CD), where B is the signal’s bandwidth, ^2 is 
the chromatic dispersion and L is the transmission distance. 
One burst consisted of one data packet and the associated 
guard interval. At the receiver, after synchronization, each 
burst is extracted and processed separately. Since the forward 
and inverse NFTs require that the signal has vanished early 
enough before it reaches the boundaries, zero padding was 
applied to enlarge the processing window as in Fig. [^c). 

IV. Simulation Results 

In this section, the performance of DBP in the NFD is 
compared with a traditional DBP algorithm based on the split- 
step Fourier method as well as with simple chromatic 
dispersion compensation; the latter works well only in the low- 
power regime where nonlinear effects are negligible. 

A. Normal dispersion 

In the normal dispersion case, a 56 Gbaud Nyquist-shaped 
transmission scheme is considered in burst mode with 256 
symbols in each packet. The duration of each packet is '-4.6ns. 
The burst size is 16000 samples and the processing window 
size for each burst after zero padding is D = 65536 samples. 
The guard time is '-10% longer than the fiber chromatic dis¬ 
persion induced memory of the link. The forward propagation 
is simulated using the split-step Fourier method (T^ with 









































































80 steps/span, i.e. 1km step. Monte-Carlo simulations were 
performed to estimate the system performance using the error 
vector magnitude (EVM) 1231 (5)]. For convenience, the EVM 
is then converted into the Q-factor 20 log(V2 erfc“^(2 BER)) 
using the bit error rate (BER) estimate (231 (13)]. 

The performance of the 100-Gb/s QPSK Nyquist-shaped 
system is depicted in Fig. as a function of the launch power 
in a 4000km link for various configurations. An exemplary 
fiber in- and output as well as the corresponding reconstructed 
input (via DBF in the NFD) are shown in Fig. It can be 
seen in Fig. that the proposed DBF in the NFD provides a 
significant performance gain of -8.6 dB, which is comparable 
with the traditional DBP employing 20 steps/span. Traditional 
DBF with 40 steps/span can be considered as ideal DBP 
in this experiment because a further increase of the number 
of steps/span did not improve the performance further. In 
the considered 4000km link, 40 steps/span DBP requires 
M = 2000 steps in total. This illustrates the advantage of 
the proposed DBP algorithm whose complexity is in contrast 
independent of the transmission distance. The performance of 
DBP in the NFD was however found to degrade rapidly when 
the launch power is sufficiently high. We believe that this 
effect can be mitigated through an extension of the processing 
window at the cost of increased computational complexity. 

A similar behavior can be observed in Fig. for a 300-Gb/s 
64QAM Nyquist-shaped system. It can be seen in Fig. that 
DBP in the NFD shows -1.5 dB better performance than DBP 
employing 20 steps/span. The difference from ideal DBP (40 
steps/span) is just --O.fi dB. This clearly shows that NFD-DBP 
can provide essentially the same performance as ideal DBP. 

B. Anomalous dispersion 

Most fibers used today (such as SSMF) have anomalous 
dispersion. Therefore, it is critical to evaluate the performance 
of the proposed DBP algorithm in optical links with anomalous 
dispersion. The DBP-NFD algorithm proposed earlier however 
cannot yet deal with signals where the function a(0, A) defined 
in 0 has roots in the upper half-plane. See Sec. |II-D[ The 
condition that the Li-norm of each burst is less than iv/2 
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Figure 5. Top: True vs reconstructed (via DBP-NFD) fiber input for a 
QPSK-Nyquist signal. Bottom: Corresponding fiber output. Only real parts 
are shown. 



is sufficient to ensure this condition. The packet duration 
in a system with anomalous dispersion thus has to be kept 
small enough in order to apply the proposed algorithm. This 
constraint is not desirable in practice, as it reduces the total 
throughput of the link because the guard interval, which is 
independent of the packet duration, must be inserted more 
frequently. The Li norm of the Fourier transform is always 
lower or equal to the Li norm of the time-domain signal. 
Motivated by this observation, OFDM was chosen instead of 
Nyquist-shaping for anomalous dispersion. 

The performances of both traditional DBP and DBP in the 
NFD are shown in Fig. DBP in the NFD can achieve a 
performance gain of -3.5 dB, which is dB better than DBP 







































Figure 7. Performance of 100-Gb/s QPSK OFDM over 4000km in fiber 
with anomalous dispersion. 



Figure 8. Li-norm and soliton-signal power ratio for lOOGb/s OFDM. 

employing 4 steps/span. The performance however degrades 
dramatically once the launch power has become larger than 
—5dBm. We attribute this phenomenon to the emergence of 
upper half-plane roots of a{0,z), which corresponds to the 
formation of solitonic components in the signal. This argument 
is supported in Fig. where the Li-norm and the ratio 
between the power in the solitonic part (E in (H p. 4320]) 
over the total signal power are plotted as functions of the signal 
power. At a launch power of —4 dBm, the signals begin to 
have solitonic components, which seems to have a significant 
impact on the DBF algorithm in the NFD. As anticipated in 
Sec. |II-D[ solitons indeed only occur above an Li-norm which 
is significantly higher than the bound 7r/2. 

V. Conclusion 

The feasibility of performing digital backpropagation in the 
nonlinear Fourier domain has been demonstrated with a new, 
fast algorithm. In simulations, this new algorithm performed 
very close to ideal digital backpropagation implemented with 
a conventional split-step Fourier method for fibers with normal 
dispersion at a much lower computational complexity. In the 
anomalous dispersion case, it was found that the algorithm 


works well only if the signal power is low enough such that 
solitonic components do not emerge. We are currently working 
to remove this limitation. 
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